Universality of spectra for interacting quantum chaotic systems 
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We analyze a model quantum dynamical system subjected to periodic interaction with an envi- 
ronment, which can describe quantum measurements. Under the condition of strong classical chaos 
and strong decoherence due to large coupling with the measurement device, the spectra of the evo- 
lution operator exhibit an universal behavior. A generic spectrum consists of a single eigenvalue 
equal to unity, which corresponds to the invariant state of the system, while all other eigenvalues 
are contained in a disk in the complex plane. Its radius depends on the number of the Kraus mea- 
surement operators, and determines the speed with which an arbitrary initial state converges to 
the unique invariant state. These spectral properties are characteristic of an ensemble of random 
quantum maps, which in turn can be described by an ensemble of real random Ginibre matrices. 
This will be proven in the limit of large dimension. 



I. INTRODUCTION 

Time evolution of an isolated quantum system can be 
described by unitary operators. Quantum dynamics cor- 
responds then to an evolution in the space of quantum 
pure states, since a given initial state \ip) is mapped into 
another pure state \tp') — U\i(j), where U = exp(— iH). 
Here H represents a Hcrmitian Hamiltonian of the sys- 
tem and the time t is set to unity. 

If the underlying classical dynamics is chaotic the 
Hamiltonian H or the evolution operator U can be mim- 
icked by ensembles of random unitary matrices [J HI ■ In 
particular, spectral properties of an evolution operator 
of a deterministic quantum chaotic system coincide with 
predictions obtained for the Dyson ensembles of random 
unitary matrices Q. The symmetry properties of the 
system determine which ensemble of matrices is appli- 
cable. For instance, if the physical system in question 
does not possess any time-reversal symmetry, one uses 
random unitary matrices of the circular unitary ensem- 
ble (CUE). If such a symmetry exists and the dimension 
of the Hilbert space is odd one uses symmetric unitary 
matrices of the circular orthogonal ensemble (COE) Q • 

If the quantum system S is not isolated, but it is cou- 
pled with an environment E, its time evolution is not uni- 
tary. One needs then to characterize the quantum state 
by a density operator p, which is Hermitian, p = p', pos- 
itive, p > and normalized, Trp = 1. Time evolution of 
such an open system can be described in terms of master 
equations Q , which imply that the dynamics takes place 
inside the set of quantum mixed states. 

The coupling of the system S with an environment E 
can heuristically be described by adding to the Hamil- 
tonian an anti-Hermitian component, H — > H' = H — 
ipWW\ where W is an operator representing the in- 
teraction between both systems The corresponding 
ensembles of non-Hcrmitian random matrices with spec- 



trum supported on the lower half of the complex plane 
were studied in @, [1| ■ For any positive value of the cou- 
pling strength parameter p the dynamics of the system 
is not unitary and eigenvalues of the evolution opera- 
tor move from the unit circle inside the unit disk - see 
Fig. [lb'. A similar situation occurs if one takes into ac- 
count dissipation in the system. Such a dynamics of 
eigenvalues of a non-unitary evolution operator in the 
complex plane was analyzed by Grobe et al. 0] and later 
reviewed by Haake 

Time evolution of an open quantum system can also be 
described in terms of a global unitary dynamics V, which 
couples together a system Sa with an other subsystem 
Sb , followed by averaging over the degrees of freedom de- 
scribing the auxiliary subsystem. Technically, the image 
of an initial state p of the system is obtained by a par- 
tial trace over the subsystem Sb, p' = Ttb[V(p ® a)V^], 
where a denotes the initial state of the environment. The 
map p' = <&(p) defined in this way is completely positive 
and preserves the trace, so it is often called a quantum 
operation [l(J, E| • Note that in this approach both inter- 
acting subsystems Sa and Sb are set on an equal footing. 
The second system, usually referred to as an 'environ- 
ment', is in fact treated symmetrically, and one may also 
consider a dual operation, in which the partial trace is 
taken over the principal subsystem Sa - compare Fig.[lfc. 

A quantum map can be described by a superoperator 
which acts on the space of density operators. If N 
denotes the size of a density matrix p, the superoperator 
is represented by a matrix $ of size iV 2 . In general such 
a matrix is not unitary, but it obeys a quantum analogue 
of the Frobenius-Perron theorem, so its spectrum is con- 
fined to the unit disk [T3|. Spectral properties of super- 
operators representing some exem plar y interacting quan- 
tum systems were analyzed in [13l4l6| . It is worth to add 
that spectra of quantum superoperators are already ex- 
perimentally accessible: Weinstein et al. [Tt} study spec- 
tra of superoperators corresponding to an NMR rcaliza- 
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FIG. 1: Schematic representation of a) an isolated quantum 
system S characterized by a Hamiltonian H and a unitary 
evolution operator U = exp(— iH); b) open quantum system 
S. The influence of an environment E can be described by 
an anti-Hermitian part of the Hamiltonian, —ifiWW*; and c) 
interacting systems S and E, in which the global evolution is 
unitary, and the non-unitarity of the evolution of Sa is due to 
the partial trace over the subsystem Sb- Panels a'), b') and 
c') show exemplary spectra of the corresponding evolution 
operators, which belong to the unit disk on the complex plane 
z = x + iy. 



of random operations [l2|- In other words, we explore the 
link between quantum chaotic dynamics and ensembles 
of random matrices. The analysis performed earlier for 
unitary quantum dynamics [l| (see Fig. QJi) is extended 
for a more general case of non-unitary time evolution of 
interacting quantum systems. This problem can be de- 
scribed by an approach closely related to the one used 
earlier to characterize quantum dissipative dynamics. To 
describe spectra of such non-unitary evolution operators 
Grobe et al. Q applied random matrices of the complex 
Ginibre ensemble |18j |. 

The key idea of this work can be visualized in Fig. [5J 
which shows a bridge established between interacting 
quantum systems, appropriate ensembles of random op- 
erations and ensembles of Ginibre matrices. Since a su- 
pcroperator describing one-step evolution operator can 
be represented as a real matrix fliij , we are goin g to ap- 
ply random matrices of the real Ginibre ensemble [20l . [2l| . 
In particular, we investigate time evolution of initially 
random pure states in a deterministic model of quantum 
baker map periodically subjected to quantum measure- 
ments, study the speed of their convergence to the invari- 
ant state and compare the results with those obtained for 
an appropriate ensemble of random operations. 




tion of exemplary quantum gates. 

For a quantum operation $ there exists an invariant 
state uj = <I>(u;). In a generic case of a typical (ran- 
dom) operation such an invariant state is unique [l2j]. 
If the action of the map is repeated n times any initial 
quantum state p converges to u> exponentially with the 
discrete time n. The rate of this convergence is governed 
by spectral properties of the superoperator, which can be 
characterized by the spectral gap, defined as the differ- 
ence between moduli of the two largest eigenvalues. 

The main aim of this work is to analyze spectra of evo- 
lution operators representing interacting quantum sys- 
tems. We demonstrate that under the condition of strong 
classical chaos and strong decoherence these spectral 
properties are universal and correspond to an ensemble 



FIG. 2: Under the assumption of strong chaos and large de- 
coherence a deterministic dynamics of a) an interacting quan- 
tum system can be described by b) an ensemble of random op- 
erations (completely positive, trace preserving maps). These, 
in turn, can be mimicked by c) random matrices of the real 
Ginibre ensemble, (which do not imply CP and TP proper- 
ties). 



This paper is organized as follows. In section II we re- 
call necessary definitions and introduce several versions 
of a deterministic quantum system model: the quantum 
baker maps subjected to measurement process. In Sec- 
tion III we introduce various ensembles of random maps. 
In section IV the spectra and spectral gaps of superoper- 
ators corresponding to baker map are analyzed and com- 
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pared with the spectra of random operations. In section 
V we continue to analyze spectra of random quantum op- 
erations and give a proof for the real Ginibre conjecture 
put forward in [l2j , and we study the fraction of the real 
eigenvalues. For large dimension N this ratio is found 
in agreement with the predictions obtained for the real 
Ginibre ensemble. 



taken over the ensemble of initially random pure states, 
po = \tp)(tp\. In the case of a generic map an exponen- 
tial convergence to equilibrium, d(t) = d(0) exp(-at) was 
reported [12j | . The convergence rate depends on spec- 
tral properties of the superoperator The spectrum 
can be characterized by the spectral gap, j = 1 — \z%\, 
which gcncrically determines the convergence rate, a = 
-ln(l- 7 ). 



II. QUANTUM OPERATIONS AND SPECTRAL 

GAP 



We shall start reviewing the necessary notions and def- 
initions. Let us define the set of quantum states Mn 
which contains all Hermitian, positive operators p of size 
N x N with trace set to unity. A quantum linear map $ 
acting on M. n is called completely positive (CP) if pos- 
itivity of the extended map, (<& <g) 1m){p) ^ : Vp holds 
for an arbitrary dimension M of the extension, and trace 
preserving (TP) if Tr($(p)) = Tr (/=>). Any CP TP map is 
called quantum operation or stochastic map. If a quan- 
tum operation $ preserves the identity, $(1/N) = 1/N, 
the map is called bistochastic. 

According to the dilation theorem of Stinespring [22| 
any CP map may be represented by a finite sum of M 
Kraus operators, 



M 



(1) 



m— 1 



If the Kraus operators A m satisfy the identity resolution, 
^ ra ( J 4" 1 ) t i m = Ijv, the map $ is trace preserving. The 
corresponding superoperator can be expressed as a sum 
of the tensor products [llj , 
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A m <g) {A m y 



(2) 



where the * denotes the complex conjugation. Let z% 
with i = 1 , . . . , N 2 denote the spectrum of $ ordered 
with respect to the moduli, \z\\ > \z2\ > ■ • • > \zn 2 \ > 0. 

A quantum stochastic map $ sends the compact, con- 
vex set A4jy of mixed density matrices into itself. Hence 
such a map has a fixed point, the invariant quantum state 
u> = $(oj). Thus the spectrum of any superoperator $ 
representing a quantum operation contains an eigenvalue 
Z\ equal to unity, while all other eigenvalues belong to the 
unit disk. In the case of unitary dynamics the leading 
eigenvalue is degenerated, but for a random stochastic 
map the invariant state is generically unique [l2j . In this 
case any pure state, (?/>|, converges to the equilibrium 
state u) if transformed several times by the map $. 

To analyze the rate of convergence to u we analyze an 
average trace distance to the invariant states, 



d(t) = (Tr\&(p )-> 



(3) 



where t denotes the discrete time (i.e. the number of 
consecutive actions of the map $), while the average is 



III. DETERMINISTIC SYSTEM: QUANTUM 
BAKER MAP SUBJECTED TO 
MEASUREMENTS 

The formalism of discrete quantum maps is applica- 
ble to describe a deterministic quantum system, period- 
ically interacting with an environment. In this work we 
concentrate on quantum dynamical systems, the classical 
analogues of which are known to be chaotic. Following 
the model of Balazs and Voros [23| we consider the uni- 
tary operator describing the one-step evolution model of 
quantum baker map, 
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? N/2 



(4) 



Here Fn denotes the Fourier matrix of size N, [Fjvljfc = 
cxp(\jk/2-KN)/\/N and it is assumed that the dimension 
N of the Hilbert space Hn is even. The standard quan- 
tum baker map B may be generalized to represent an 
asymmetric classical map, 



B K = F 



N 



F, 



N/K 







Fn-n/k 



(5) 



where K > 2 is an integer asymmetry parameter chosen 
in such a way that the ratio N/K is integer. The stan- 
dard model, obtained in the case K = 2, corresponds to 
the classically chaotic dynamics characterized by the dy- 
namical entropy H equal to In 2 [24|. This system can 
be considered as a 2-dimensional lift of an 1-dimensional 
non-symmetric shift map, 

Kx/{K - 1) :x£[0, (K - l)/K] 
f K {x) = { . (6) 

Kx-K + 1 : x G ({K - 1) /K, 1] 

Chaos in such a system can be characterized by its dy- 
namical entropy, equal to the mean Lyapunov exponent, 
averaged with respect to the invariant measure of the 
classical system. Since the uniform measure is invariant 
with respect to the map fx, the dynamical entropy h is 
equal to the mean logarithm of the slope dfx/dx, 



h{K) 



llnA-+^ 
K K 



V K 



K - 1 



(7) 



The entropy is maximal in the case K = 2, while in 
the limit K — > 00 the entropy tends to zero. Hence the 



4 



larger value of the parameter K is, the weaker chaos in 
the classical system becomes. 

In the case of the quantum system acting on the TV- 
dimensional Hilbert space the largest possible value of the 
asymmetry parameter reads K = N . Thus the limiting 
case of the classically regular system cannot be obtained 
for any finite N. The limit of vanishing dynamical en- 
tropy, h — >• 0, can be approached only in the classical 
limit TV — > oo of the quantum system. 

A generalized variant of a non-unitary baker map in- 
troduced by Saraceno and Vallejos described a dissipa- 
tive quantum system [25j . In this work we study another 
model of non-invertible quantum baker map analyzed in 
[HIHH, which is deterministic, conserves the probability, 
and is capable to describe projective measurements or a 
coupling with an external subsystem. 

In general there exist M different outcomes of the mea- 
surement process and thus the map is described by a col- 
lection of M Kraus operators. The simplest nontrivial 
case of M = 2 corresponds to dividing of the phase space 
into two parts, which we can choose to be the 'lower' and 
the 'upper' part. Such a measurement scheme allows one 
to write down the quantum operation corresponding to 
the 'sloppy baker map', in which both pieces of the clas- 
sical phase space are not placed precisely one by another, 
but in each step an overlap of a positive width takes place. 
In the classical model the upper piece of the phase space 
is shifted down by A/2 - see Fig. [3J; - so the invariant 
measure lives in the rectangle of the width (1 — A). 

1/K 
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FIG. 3: Sketch of the classical dynamical system acting on 
the torus: b) reversible asymmetric (K — 4) baker map Bk, 
c) irreversible sloppy baker map Bk,a in which in each step 
the upper part of the phase space is shifted down by A/2, d) 
double sloppy baker map Bk,a,a in which both parts of the 
phase space are shifted vertically by A/4. 

To represent the shift in the quantum analogue of the 
map one uses a unitary translation operator V such that 
V N = Ijv and any momentum cigenstate \k) is shifted 
by one, V\k) = \k + 1). Hence the shift down by A/2 
is realized by the unitary operator, V~ NA / 2 . Thus the 



stochastic map describing the quantum sloppy map [13| 
®b k ,a (p) = D b B K pB* K Dl + D t B KP B K D\. (8) 

consists of two Kraus operators, which act on the unitar- 
ily rotated state p' = BkpB^ k , 
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D t = V- NA / 2 Fl 
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(9) 



The operator D b describes the projection on the lower 
part of the phase space, while the definition of the oper- 
ator D t includes also the operator representing the shift 
of the upper domain down by A/2. Observe that the pa- 
rameter A may take any real value from the unit interval 
[0, 1]. However, the case A = corresponds to the baker 
map without the shift but with a measurement, so it does 
not reduce to the standard unitary baker map Bk- 

One can also consider another classical model of double 
sloppy map, in which both domains are simultaneously 
shifted by A/4 towards the center of the phase space [27| 
- Fig. [3ji. To write down the corresponding quantum 
model Br a.a one needs thus to modify both Kraus op- 
erators, D b -U V NA l A D b and D t -> V~ NA ' A D t . 

Both variants of the model can be further generalized 
by allowing for a larger number M of measurement oper- 
ators, represented by projectors on mutually orthogonal 
subspaces. For simplicity we assume here that the di- 
mensionalities of all these subspaces are equal and read 
N/M. Varying the parameter M one may thus control 
the degree of the interaction of the baker system with the 
environment and study the relation between the decoher- 
ence in the interacting quantum system and the spectrum 
of the corresponding superoperator. 

Increasing the asymmetry parameter K one can de- 
crease the degree of classical chaos. To increase the de- 
gree of chaos one may just apply the quantum baker 
map twice, since the classical dynamical entropy of such 
a composite map is equal to 2 In 2. In general one can 
allow for an arbitrary number of L of unitary evolu- 
tions, and replace unitary B by B L . Alternatively one 
can say that the non-unitary measurement operation is 
performed only once every L periods of the unitary evo- 
lution. Choosing the parameter L to be of order of N 
one can assure that the quantum dynamics is as 'chaotic' 
as allowed by the quantum theory, what can be quantita- 
tively characterized by the quantum dynamical entropy 
[11433. 

Therefore the generalized model of quantum sloppy 
baker map we are going to analyze here depends on the 
classical asymmetry parameter K, the width of the clas- 
sical shift A, the number of free evolutions L, and the 
quantum parameter M denoting the number of measure- 
ment operators, 



4> 



,{ P ) = J2 d ^(b k ) l p(bI) l dI. (10) 
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Additionally, for each set of parameters of the model one 
can choose the appropriate set of projection operators 
D m which correspond to the shift applied on one or on 
two parts of the classical phase space. 

Note that the measurement process can also be inter- 
preted as an interaction with a measurement apparatus, 
described by an auxiliary Hilbcrt space of M dimensions. 
Thus the model (|10[) represents an interacting quantum 
system and belongs to the general class of quantum oper- 
ations defined by |T]). A rich structure of the model and 
the possibility to tune independently several parameters 
of the quantum system allows us to treat this model as a 
valuable playground to investigate spectral properties of 
supcroperators, which represent non-unitary dynamics of 
interacting quantum systems. 
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The spectra of the superoperator of the generalized 
sloppy baker map (fTU|) were found to depend weakly on 
the shift parameter A. However, other parameters of the 
model (namely N, K, L and M) influence properties of 
the spectrum considerably - see Fig. [4] 

As the asymmetry parameter K increases the differ- 
ence between the sizes of two domains which form the 
classical phase space becomes larger. In the extreme 
limit of K — > N — >• oo the classical system becomes only 
marginally chaotic, the eigenvalues are attracted to the 
unit circle and the spectral gap 7 = 1 — \z%\ disappears. 

On the other hand, if we increase the degree of the clas- 
sical chaos by increasing the number L of unperturbed 
unitary evolutions, the size of the spectral gap does not 
change, but the spectrum fills the complex disc of radius 
R = I Z2 1 almost uniformly. Eventually, an increase of the 
number M of the measurements results in a faster deco- 
herence in the system. This is reflected by an increase 
of the spectral gap 7. In fact the radius R = 1 — 7 of 
the disk supporting the spectrum decreases with M as 
M~ x l' 2 . This observation - demonstrated in Fig. [7]- will 
be explained in section IVl 

IV. ENSEMBLES OF RANDOM OPERATIONS 
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FIG. 4: Exemplary spectra of the evolution operator for the 
sloppy baker map for several values of the parameters of the 
model. The dimension of the Hilbert space N = 64, param- 
eter M = 2, and the shift width A = 1/4 are kept fixed. 
A generic spectrum for K — 4, L = 16 is shown on panel 
b). The subplots a) and c) are obtained for the cases of a 
weak classical chaos for K — 64, L = 1 and K = L = 32 re- 
spectively, while the last case d) shows the spectrum for the 
double sloppy map Bk,a,a for K = 64 and L — 64. 

We constructed quantum operations representing the 
generalized sloppy baker map (fT0|) for several sets of the 
parameters of the model. In each case the superopera- 
tor $ was obtained according to the expression ^ and 
diagonalized to yield the complex spectrum belonging to 
the unit disk. 

In the case of several measurement operators, M > 2, 
the quantum baker map represents a non- unitary dynam- 
ics. Under the condition of classical chaos the leading 
eigenvalue z\ = 1 is not degenerated and all remaining 
eigenvalues are located inside the disk of the radius equal 
to the modulus of the sublcading eigenvalue R — \z2\. 



In this section we propose three different ensembles 
of random stochastic maps acting on the space M jy of 
mixed states of size N with different physical interpreta- 
tion. We assume that all unitary matrices U used below 
are drawn according to the Haar measure on the unitary 
group of corresponding dimension unless stated other- 
wise. 

1. Environmental representation of a random stochas- 
tic map [l2l | . Choose a random unitary matrix U of com- 
posite dimension NM and construct a random map as 

^E{p) = r rr E [U{p»\v){v\)U^]. (11) 

It is assumed here that the environment, initially in an 
arbitrary pure state \v) is coupled with the system p by a 
random global unitary evolution U. The stochastic map 
is obtained by performing the partial trace over the M- 
dimensional environment. 

2. Random external Gelds defined as a convex combi- 
nation of M unitary evolutions Q 

M 

Mp) = p™ u ™p u ™ ( 12 ) 

m— 1 

where p m are positive components of an arbitrary proba- 
bilistic vector of size M, Y^m=iP m = ^ ^ unitaries 
U m G U(A) are independent random Haar matrices. 
Random external fields form an example of bistochas- 
tic maps. They represent the physically relevant case in 
which the quantum system is subjected randomly with 
one of M given unitary operations and can also be inter- 
preted as quantum iterated function systems |3lj . 



6 



3. Projected Unitary Matrices acting on states of a 
composite dimension, TV = KM. All M Kraus opera- 
tors are formed by unitarily rotated projection operators, 
A m = P m U for m = 1, . . . , M which leads to the map 

M 

<S> p (p) = P m U P U^P m , (13) 

m— 1 

where U is a fixed random unitary matrix. Here P m = 
PJj = P 2 n denote projective operators on K dimensional 
mutually orthogonal sub-spaces, which satisfy the identity 
resolution, ® P m = Ijy. This ensemble of bistochastic 
maps corresponds to a model of deterministic quantum 
systems, in which unitary dynamics is followed by a pro- 
jective measurement. 
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FIG. 5: Spectra of superoperators corresponding to typical 
random maps generated according to a) the ensemble $e, and 
b) ensemble $p (for definitions see Sec. II). All maps act on 
quantum states of size N = 64, while the parameter of the 
model reads M = 2. Note that in all cases the spectrum is 
contained in the disk of radius R = l/y/2 apart of the leading 
eigenvalue marked by 'x'. 



it from unity (unitary dynamics) to TV 2 , which describes 
a generic random stochastic map. 

We have generated several realizations of random maps 
from the ensembles $e and <&p introduced above. Ex- 
emplary spectra of superoperators for maps pertaining 
ensembles obtained for M = 2 arc shown in Fig. [5] In 
the latter case we superimposed the spectra from two 
realizations of the map $p since by construction N 2 /M 
eigenvalues of the superoperator are equal to zero. 

In general, the spectra of random maps could be used 
to describe the spectra of deterministic system (|10p under 
the condition of classical chaos and large decoherence. 

Numerical results performed for various models of 
quantum maps reveal an exponential decay of the mean 
trace distance ([3]) to the invariant state. A comparison 
of such a time dependence of the mean trace distance d 
for quantum baker maps and random quantum maps is 
shown in Fig. [5] Interestingly, for a fixed system of size N 
the convergence rate a increases with the number of the 
measurements as a ~ ^ In M , but it seems not to depend 
on the assumption, whether the global random evolution 
matrix U is taken from the orthogonal or unitary circular 
ensemble. 




In the ensembles of random maps defined above the 
integer number M 1 determines the number of Kraus 
operators and serves as the only parameter of each en- 
semble of random maps. Observe that in the special case 
M = 1 the dynamics reduces to the unitary evolution, 
so both variants of the model are used to describe quan- 
tum systems with or without a generalized antiunitary 
symmetry [l|. 

As shown in QjJ the flat measure in the set of stochas- 
tic matrices is obtained for the coupling of the system 
with an environment of dimension M = N 2 , so that 
the Choi matrix, 73$ := ($ ® t)\cf> + ) (4> + \ of size N 2 

has full rank. Here \4>+) = jj Y2%=i I*)*) represents the 
maximally entangled state on the bipartite Hilbert space, 
T~Ln <8> T-Ln- Due to the theorem of Choi the condition of 
complete positivity of the map is equivalent to positivity 
of the Choi matrix, 

$ is CP & 79$ > . (14) 

In general the discrete parameter M characterizes the 
strength of the non-unitary interaction and we shall vary 



FIG. 6: Time dependence of the average trace distance d of 
random pure states to the invariant state for a) exemplary 
random quantum maps $e and b) quantum baker map with 
parameters: K = 4, L = 16, A = 1/4. Figure drawn in a log- 
scale shows the exponential decay for dimension N = 24 and 
parameters M = 8(o) and M = 12 (a). The average is taken 
over the set of 16 initial projectors. 

Further numerical investigations confirm an expected 
relation between the size rate of the convergence of a typi- 
cal quantum state towards the fixed point of the map and 
the spectral gap. As shown in Fig.[7]obtained for random 
operations as well as the generalized quantum baker map 
the radius R = 1 — 7 of the disk in the complex plane, 
which contains all but the leading eigenvalue, decreases 
with the number of measurements as 

R =H~J=. (15) 

In the next section we present an explanation of this re- 
lation based on the theory of random matrices. 
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FIG. 7: Modulus of the subleading eigenvalue R — \z2\ as a 
function of the dimension M of the environment for random 
quantum operation $e for TV = 4 ( x ) and for the sloppy baker 
map for N = 64, K = 4, L = 16, M = 2, A = 1/4 (■). Solid 
line shows the fit according to eq. (|15[) . 



V. QUANTUM OPERATIONS AND REAL 
GINIBRE ENSEMBLE 



Looking at the spectrum of a random stochastic map 
$e shown in Fig. [5] one can divide the entire spectrum 
into three parts: i) a single eigenvalue z\ = 1; ii) Af n real 
eigenvalues distributed along the real line with a density 
P^(x), hi) remaining A/" c complex eigenvalues Zi, the 
distribution of which can be described by a density P c (z) 
on a complex plane. 

Any density operator p of size N can be represented 
using the generalized Bloch vector representation 



AT -1 

£■ 

i=0 



(16) 



Here A* denotes the generators of SU(7V) such that 

p\ a l G K for 
a N 2_ 1 ] is called 
the generalized Bloch vector. Thus 



Tr(X i X>) = S ij and A cx 1. Since p 
i = 0, . . . , N 2 — l. The real vector [ao, • 



(17) 



The Bloch vector can also be used to represent an ar- 
bitrary operation. Using Kraus operators A m one repre- 
sents the element <fr lJ of the superoperator $ in a form 

$ y = Tr (A 1 $ (X j )) = Tr ^ A 1 A m A 3 (A m ) f , (18) 



where i,j = 1, . . . , -ZV 2 — 1. This square matrix of order 
TV 2 — 1 will be called C. In a similar way we introduce 



the vector k and hnd that the remaining elements of the 
matrix $ do vanish, 

$ io = Tr^ AM m A° {A m ) f = Tr A* A ^ a m {A m ) f = (n) t 

ra rn 

(19) 

m 

(21) 

Hence the superoperator $ can be represented as a real 
asymmetric matrix 



1 

K C 



(22) 



where the N' := N 2 — 1 dimensional vector k represents 
a translation vector while the N' x N' real matrix C is 
a real contraction 
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Thus the spectrum of $ consists 
of the leading eigenvalue, equal to unity, and the spec- 
trum of C. Note that the complex eigenvalues of the real 
matrix C appear in conjugated pairs, z and z, which is 
a consequence of the fact, that the map $ sends the set 
of Hermitian operators into itself, so as seen above the 
superoperator can be represented by a real matrix fl9l ]. 
Since a map acting on states of size N is represented by 
a superoperator of dimension Af 2 the following normal- 
ization relation holds, 1 + 7V R + 7V C = iV 2 . 

In the case that C has only real eigenvalues one can 
bring C by an orthogonal transformation O to lower tri- 
angular form 



c = o{~+ A)er 



(23) 



where A = diag(zi, . . . , zn>) while 3 has elements only 
below the diagonal. Thus 



dC = Op^dOiZ 



A) - (H + A)Q- 1 dO- 
+ dZ + dA}0- 1 . 



Hence the measure DC is given by 



DC = 



j( Zi ^ Zj ) n d Zk Hio^do)^ n <e, 



i<j 



i<j 



3<i 



(24) 



(25) 



where the Vandcrmonde determinant is the Jacobian 
of the transformation from (0~ l dO)ij Aj to (0~ 1 dO)ij. 
Thus the measure dp(A) has the form 



Dp(X) = 



i<j 



e(i?$ ^ o) 



(26) 



where in the last factor the positivity conditions of the 
corresponding Choi matrix is averaged by integrated over 
the measure DkDE ]JO~ 1 dO. This factor is expected to 
be a smooth function of the eigenvalues z\, . . . , znl In 
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the case that C has a certain number of complex conju- 
gate eigenvalues D/i(A) is of similar form, but the prod- 
uct of differentials dz k has to be interpreted as exterior 
product [2l|. It turns out that for large dimension N' 
the measure d/j,(A) is given by the real Ginibre ensemble 
with the bulk of eigenvalues inside a certain disk in the 
complex plane. To prove this let us go back to the ma- 
trix representation of $ in terms of M Kraus operators 
A m ,m = 1,...,M. Then 

M 



m—1 



kl 



where A^t are the matrix elements of Kraus opera- 
tors and * denotes the complex conjugation; i,j,k,l = 
1, . . . , N. The Kraus operators obey 



M N 



ki 



(28) 



m—1 i—1 



thus it is natural to assume that A™, represent N columns 
of a matrix U drawn from a circular unitary ensemble of 
dimension NM i.e. U € D(NM). Using formulas by 
Mello [13 for the first four moments of U(iVM) we are 
able to find exactly the first two moments of matrix el- 
ements For example for U € U(JV) : (UbpU* a ) = 
bab&up/N. Here (• • • ) means the averge over the unitary 
group. This implies here: 



1 



M 

($ ij!kl ) = J2(Ar k (A? l Y) = jj6 ji S lk . 

771—1 



(29) 



In this way we can derive the exact second moments: 

(NM) 2 - 1 ^ M2 ^ Sr J Skl5 ki + MS ii s fj S kk S u) 

~ NM (NM) 2 -1 ( M2Si ^ 5 kk 5 n + MS u 5 fMhi) ■ 

(30) 

We see that in the limit of large N the first two cumulants 
are identical to those of the Gaussian distribution (with 
variance 1/(N 2 M)) 



P($) cx cxp NM *i 



N 2 M 



,kk 



ijkl 



(31) 

The factor 1/2 is due to the symmetry property 3>y,fcj = 
^jilk- We will argue below that for large M in addi- 
tion the higher cumulants can be neglected. Hence the 
superoperator $ associated with a random map can be 
described (up to the one eigenvalue 1) by the real Ginibre 
ensemble with eigenvalues inside a disk of radius 1 / VM, 
where M is the number of random Kraus operators defin- 
ing $. This can also be seen by going back to the real 
matrix representation (TTJ 



Let us now argue that for large M we can neglect 
higher cumulants. First of all for large TV the elements 
A"^, forming a minor of U € U(iVM), are essentially inde- 
pendent Gaussian variables with zero mean and variance 
1/NM. Thus as a consequence of the central limit theo- 
rem for large M &ij,ki as sum of M essentially indepen- 
dent identically distributed variables is again Gaussian 
with variance M/(NM) 2 = l/N 2 M. Also in the bulk 
the different matrix elements of $ arc independent. The 
average of is given by 5 i:j 8 kl /N. 

To investigate the density of complex eigenvalues z = 
x + iy of the superoperator $ we analyzed their radial 
probability distribution P(r), where r = \z\. The real 
eigenvalues are taken into account for this statistics. Fig. 
[5] shows a comparison of numerical data obtained for sev- 
eral realizations of quantum baker map, projective ran- 
dom operations and real random matrices pertaining to 
the Ginibre ensemble. The data are represented in the 
rescaled variable tm = ry/M so that the radius of the 
disk of eigenvalues is set to unity. In all three cases 
displayed in the figure the radial density grows linearly, 
which corresponds to the flat distribution of eigenvalues 
inside the complex disk, in agreement with the predic- 
tions of the Ginibre ensemble. These results obtained for 
A^ = 32 show a smooth transition of the density in the 
vicinity of the boundary of the disk at tm = 1, which be- 
comes more abrupt for larger N. In the asymptotic case 
A^ — > oo, the density of rescaled eigenvalues is described 
by the circular law of Girko, 



p c (z)~e(i-|z| 



(32) 



derived for complex Ginibre matrices. 

The spectra of real random Ginibre matrices display 
a more subtle structure. A finite fraction of all eigen- 
values are real, in analogy to the mean number of real 
roots of a real polynomial [3^, [34[ . Real eigenvalues of a 
real Ginibre matrix cover the real axis with a constant 
density. Furthermore, for large dimension the density of 
complex eigenvalues is known to be asymptotically con- 
stant in the complex disk except for a small region near 
the real axis (2pj . 

To check for what random operations these effects can 
be observed in the spectrum of the superoperator we an- 
alyzed the average number (A/' IR )<i> of real eigenvalues of 
the superoperator $. For any realization of $ we have 
Af R = #REAL/(A^ 2 - 1) where #REAL is the number of 
real eigenvalues of the real matrix of size N' = N 2 — 1 . 
These data are compared with predictions for the real 
Ginibre ensemble, hereafter denoted by (A/" r )rg- The 
following expression for the mean number of real eigen- 
values of a real Ginibre matrix of size A^ 2 — 1 was derived 
in (35rS3 



R\ 



IRG 



= 1 



V2 [ 1 t 1 ' 2 (l-t N '- 2 )dt 



7T J Q (l-t)3/2( 1+f ) 



-Vn 2 

IT 



1 as N 



(33) 
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FIG. 8: Radial density of complex eigenvalues of spectra of su- 
peroperators corresponding to a deterministic model of sloppy 
baker map (x), projective random operations (o) and the 
spectra of real random matrices of the Ginibre ensemble (■). 
The size of each matrix is TV 2 = 32 2 , the number of Kraus op- 
erators is M = 16 so the density is shown as a function of the 
rescaled radius tm '■= r^/M = 4r. The tail of the distribution 
beyond the point tm — 1 (equivalent to r = 1/4) does not vi- 
olate therefore the quantum analogue of the Frobenius-Perron 
theorem. 



These analytical results suggest to introduce a rescaled 
ratio 



■q := 



Viv 2 - 1 



(34) 



to make easier a comparison of data obtained for vari- 
ous systems of size N. Numerical results presented in 
Fig. [S]show that the superoperators associated with ran- 
dom maps are characterized by a non-zero fraction of 
real eigenvalues. In the case of strong interaction with 
the ancilla of the size M = N 2 the dynamical matrix D<j> 
has full rank and the rescaled fraction of real eigenvalues 
of $ coincides with the prediction for the real Ginibre 
ensemble. 

To demonstrate further spectral features characteris- 
tic of the real Ginibre ensemble we analyzed spectra of 
superoperators and investigated the cross-section of the 
probability distribution P(z) along the imaginary axis. 
Numerical data of this distribution denoted as P*p (y) ob- 
tained for an ensemble of random maps $e acting on the 
states of size N = 3 are shown in Fig. [TO] 

In order to compare these data with predictions of 
the real Ginibre ensemble we need to assure a suitable 
normalization. Let us rescale the imaginary axis as 
y <— > yu = V~My, so that the rescaled formula (25) of 
1371 takes the form 



2M 



cxp (2My 2 ) Merfc(|y| 



(35) 



1 

0.95 
0.9 

0.85 
0.8 
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FIG. 9: Rescaled ratio of the real eigenvalues r\ of the superop- 
erator for random operations with M = iV 2 (□), M = N (o) 
and real Ginibre matrices (a) as a function of the matrix size 
N. Solid horizontal line at W2/tt represents the asymptotic 
value of the normalized ratio implied by (|33[) . 



Making use of the standard estimations 



1 



< exp(x 2 ) / exp(-t 2 )dt < 



+ \^ 2 + 2 

(36) 

(see formula (7.1.13) at page 298 of Abramowitz and Ste- 
gun HU), and the definition of the complementary error 
function erfc(z), one obtains from eq. (|35|) an explicit 
form for lower and upper bounds for the rescaled distri- 
bution in the vicinity of the real axis 



1 



1 



1 



l 



1 

^ -- 



1 



1 + 4 



2 1 
7T My' 2 



(37) 

As shown in Fig. [10] these bounds arc rather precise and 
describe well the numerically observed density Pf'(y) 
of complex eigenvalues of the superoperators along the 
imaginary axis. 



VI. CONCLUDING REMARKS 

In this work we analyzed spectra of non-unitary evo- 
lution operators describing exemplary quantum chaotic 
systems and the time evolution of initially pure states. 
We have chosen to work with a generalized model of 
quantum baker map subjected to measurements [l3l . [27| , 
which allows one to control the degree of classical chaos 
and the strength of the interaction with the environment. 
The size of the quantum effects, proportional to the ratio 
of the Planck constant to the typical action in the sys- 
tem, is controlled by the size N of the Hilbcrt space used 
to describe the quantum system. The classical limit of 
the quantum model corresponds to the limit N — > oo. 
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0.35 
P c (y ) 

1 Vy M' 
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0.15 

0.05 





FIG. 10: Numerical data for the density -Pi C (j/a/) of complex 
eigenvalues along the rescaled imaginary axis obtained for an 
ensemble of random operations <3>e of dimension N = 8 and 
with parameter M = N 2 = 64 (■) are compared with lower 
and upper bounds (|37p obtained for small \y\ from the real 
Ginibre ensemble and represented by thick lines. 



Due to a quantum analogue of the Frobenius-Perron 
theorem the evolution operator has at least one eigen- 
value equal to unity, while all other eigenvalues are con- 
tained in the unit disk in the complex plane. In a generic 
case the leading eigenvalue is not degenerated and the 
corresponding eigenstate represents the unique quantum 
state invariant with respect to the evolution operator. 

Investigating the time evolution of initially random 
pure states we found out that in a generic case they con- 
verge exponentially fast to the invariant state. The rate 
of this relaxation to the equilibrium depends on the size 
of the spectral gap, equal to the difference between the 
moduli of the first and the second eigenvalues of the evo- 
lution operator. In particular, the relaxation rate a was 
found to depend on the number of measurement opera- 
tors M as | In A/. 

Spectral properties of evolution operators of determin- 
istic quantum systems interacting with the environment 



were compared with spectra of suitably defined ensembles 
of random matrices. Note that an idea to apply random 
matrices to model evolution operators of open determin- 
istic quantum systems was put forward by Peplowski and 
Haake [39j , but random maps used therein are not neces- 
sarily completely positive. This property is by construc- 
tion fulfilled by the ensemble of random stochastic maps 
introduced in [12( and by two other ensembles of random 



bistochastic maps used in this work. 

In general, the spectra of non-unitary operators cor- 
responding to quantum deterministic systems display a 
wide variety of structures, which depend on classical pa- 
rameters as the degree of chaos of the corresponding clas- 
sical system characterized quantitatively e.g. by its dy- 
namical entropy. The spectra depend also on quantum 
parameters as the dimension of the Hilbcrt space and the 
character of the interaction with the environment, which 
governs the strength of the decoherence effects. However, 
under an assumption of strong classical chaos and a uni- 
form coupling of the system analyzed with all the states 
of the M-dimensional environment the spectral proper- 
ties the corresponding evolution become universal: The 
spectrum consists of a leading eigenvalue equal to unity, 
while all other eigenvalues cover the complex disk of ra- 
dius R = 1/VM. 

In the asymptotic limit N — > oo the density of complex 
eigenvalues becomes uniform in the disk, besides the re- 
gion close to the real axis. As the size of the environment 
M is equal to N 2 , which implies strong decoherence, the 
dynamical matrix D$ describing the quantum map $ has 
full rank, so it can be considered as generic. In this very 
case the spectral statistics of this region of the complex 
spectrum of the superoperator and the fraction of its real 
eigenvalues coincides with predictions of the reai Ginibre 
ensemble, a proof of which is given in this work. 
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